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The study of the production of heavy quarkonium is crucial for a thor- 
ough understanding of Quantum Chromodynamics (QCD). This note re- 
ports the measurements of the J/ip, Xc and double charm production cross 
section, and discusses the prospects for the J/ip polarization at LHCb. 
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is 1 Introduction 



14 Understanding the charmonium hadroproduction mechanism has been a long-term 

15 program both experimentally and theoretically. At the LHC two different compo- 

16 nents contribute: the "prompt" component, which includes "direct" production in 

17 the pp collisions and the "feed down" from higher charmonium states, and the "de- 
is layed" component, coming from the 6-hadron decays. In the direct production the 

19 cc pairs are expected to be created mainly through gluon-gluon fusion at the lead- 

20 ing order (LO) and the bound states are formed in the final color-singlet states. The 

21 most recent models allow the formation of the cc pairs also through color-octet states, 

22 which evolve toward the final state via exchange of soft gluons. Such evolution is de- 

23 scribed in terms of a non-perturbative QCD (NRQCD) factorization approach. 

24 The LHCb collaboration has given many contributions in understanding the quarko- 

25 nium production mechanism, some of them reported here. The measurement of the 

26 J I ip production cross section together with the status of the polarization analysis will 

27 be presented. Both measurements provide a critical test for the color-singlet [TJ |2] 
2B and color-octet models [3]. The study of Xc, double J /if) and J/ip production asso- 

29 ciated with an open charm hadron will also be presented. In particular the last two 

30 are rare processes and they can be useful to investigate the contributions from other 

31 mechanisms, such as the Double Parton Scattering (DPS) [I1IB1IH]. 

32 2 The LHCb detector 

33 The LHCb detector [10] is a single-arm forward spectrometer covering the pseudo- 

34 rapidity range 2 < r\ < 5, designed for the study of particles containing b or c 

35 quarks. The detector includes a high precision tracking system consisting of a silicon- 

36 strip vertex detector surrounding the proton-proton interaction region, a large-area 

37 silicon-strip detector located upstream of a dipole magnet with a bending power of 

38 about 4 Tm, and three stations of silicon-strip detectors and straw drift-tubes placed 

39 downstream. The combined tracking system has a momentum resolution Sp/p that 

40 varies from 0.4% at 5 GeV/c 2 to 0.6% at 100 GeV/c 2 , and an impact parameter 

41 resolution of 20 /^m for tracks with high transverse momentum. Charged hadrons 

42 are identified using two ring-imaging Cherenkov detectors. Photon, electron and 

43 hadron candidates are identified by a calorimeter system consisting of scintillating-pad 

44 and pre-shower detectors, an electromagnetic calorimeter and a hadron calorimeter. 

45 Muons are identified by a muon system composed of alternating layers of iron and 

46 multiwire proportional chambers. The trigger consists of a hardware stage, based 

47 on information from the calorimeter and muon systems, followed by a software stage 
4B which applies a full event reconstruction. 
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49 3 J/ip cross section measurement 



50 The cross section is measured selecting J/ip decaying to two muons: the data sample 

51 corresponds to an integrated luminosity C = (5.2 ± 0.5) pb _1 of pp collisions at y/s = 

52 7 TeV recorded by the experiment during September 2010. The double differential 

53 cross section in J / ip px and y is defined as following: 

(fa N(J/ip -> 



dp T dy L x e tot x B(J/ip ->■ fi + n~)Ap T Ay 

54 where N(J/ip — > yU + /i~) is the number of selected J/ip decaying in two muons, 

55 C is the integrated luminosity, etot is the total efficiency (estimated from Monte 

56 Carlo including the detector acceptance, the reconstruction and trigger efficiency), 

57 B(J/ip —> /i + /i~) is the branching ratio of the J/ip — > decay, Apt and Ay 
5B are respectively the J/ip transverse momentum and rapidity bin sizes. The analysis 
59 selection requires at least one reconstructed primary vertex in each event. The J/ip 
eo candidates are formed from pairs of opposite sign charged tracks reconstructed in the 
ei tracking system and identified as muons. The two muons must have a good quality of 

62 the track fit and originate from a common vertex. To separate the prompt and the de- 

63 layed component the J/ip pseudo proper time is used, defined as t z = <yZj/,p z ^ mj i^ ; 

64 where zj/$ and zpy are the J/ip decay vertex and the primary vertex positions along 

65 the beam axis and mj/^ and p z are respectively the mass and the momentum com- 
ee ponent of the J/ip along the beam axis. 

67 Figures [U and [2] show respectively the double differential prompt cross section 

es and the differential prompt cross section integrated over rapidity as a function of 

69 pt- Results are compared with the prediction of three different theoretical models 

70 (Colour Singlet Model, Colour Octet Model and Colour Evaporation Model). In FigfS] 

71 the double differential cross section of the delayed component is shown and in FigH] 

72 it is integrated over rapidity and compared with the FONLL computation. The total 

73 integrated cross sections are 

^prompt = (10.52 ± 0M(stat.) ± lA0(sys.)±H*(pol))iib, (2) 

Vfromb = [1.14±0.01(sfat.) ±0.16(sys.)}iib. (3) 

75 The first and the second uncertainties are the statistical and the systematic, where 

76 the main sources of systematic uncertainty come from the luminosity measurement, 

77 the tracking and trigger efficiency. The third uncertainty on the prompt cross section 

78 is due to the unknown polarization of the J / ip and it is estimated calculating the total 

79 efficiency in two possible extreme scenarios of fully transverse and fully longitudinal 
so polarization. The deviation from the case of zero polarization is assigned as systematic 
si uncertainty to the measurement. From Eqj3] the bb cross section is extrapolated to 
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Figure 3: Double differential cross sec- 
tion of the delayed J ftp component. 



Figure 4: Differential cross section of 
delayed J/ip component compared with 
FONLL computation. 
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82 the full solid angle using the formula 



a{pp — > bbX) = a 47r 



2B{b J/ipX) 



[288 ± 4(stat.) ± A8(sys.)} fib. 



(4) 



83 All these results have been published in Ref . [TT] . 

84 4 Outlook for polarization measurement 

85 With the full 2011 data sample the prompt J/ip polarization will be measured study- 
se ing the full angular distribution of the two muons: 



87 where 9 and <fi are the polar and azimuthal angles in the helicity frame (using the 

as J/ip momentum as polarization axis). The measurement will be performed in bins of 

89 J /ip transverse momentum and rapidity. The statistical sensitivity should be under 

90 0.15 for Xg and about 0.01 for A^ and Xg^. The systematic uncertainty is expected to 

91 be of the same order of magnitude. 

92 5 Xc production 

93 The study of the J/ip production through the radiative decays of the \ c states pro- 

94 vides a useful test of both the color-singlet and color-octet model. Moreover it is 

95 fundamental for the J/ip polarization measurement, as the directly produced J/ip 

96 and those coming from \c decays can carry different polarization and this represents 

97 a possible source of uncertainty for the polarization measurement of the prompt com- 

98 ponent. The measurement of the fraction of J/ip coming from Xc decays can quantify 

99 this uncertainty. 

100 The relative cross section cr(x c2 )/o"(x c i) is measured using two different data sam- 

101 pie acquired by the LHCb experiment during the 2010 and 2011, respectively of 37 

102 pb _1 [12] and 370 pb" 1 [13]. In both cases the \ c states are identified through their 

103 radiative decay Xc J 7V7 with the J/ip decaying to two muons J/ip -> /x 4 \T . For 

104 the first measurement, made with a smaller data sample, the photons reconstructed 

105 in the calorimeter system have been used. This allows to have a higher statistics but 

106 the poor resolution of the calorimeter doesn't permit to separate the two xa and Xc2 

107 states. In the second measurement the photons converted in the detector material 
us before the magnet have been used, 7 —> e + e~. In this way it is possible to take ad- 
109 vantage of the good resolution of the tracker, which allows to resolve the two states, 
no as it is shown in FigJSJ 



dN 



oc 1 + Xg cos 2 9 + X^ sin 2 9 cos 20 + sin 29 cos </>, 



(5) 
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Figure 5: Invariant mass difference 2 * 6 8 10 1 p t (j/ v > [Gev/c] 
spectrum Am = m(x c ) — m(J/if)) = 



m(e+e fi+fi ) - m(/i + ii ) using con- Figure 6: Relative cross section 
verted photons. <T(x<fl)MXci) in bins of J /iff trans- 

verse momentum. 

in In both measurements the efficiency is determined from the Monte Carlo simulation 

n2 and the number of signal events is extracted with a fit to the invariant mass difference 

in spectra, in four bins of J /if? transverse momentum. 

in The results for the relative cross section u(Xc2)/o'(Xci) are shown in FigfBJ In the 

us plot the inner error bars correspond to the statistical uncertainties and the outer bars 

lie correspond to the sum of all the sources of systematic uncertainties. The shaded 

ii7 area represents the maximum effect due to the unknown \ c polarization. In green 

us and black the results obtained reconstructing the photons in the calorimeter (2010 

n9 statistics) and the converted photons (2011 statistics) are shown. The results from 

120 the CDF collaboration are shown in magenta [14] . The blue and red shaded area 

121 correspond to the color-singlet and NRQCD prediction respectively. 

122 The Xc to J /iff ratio has been measured with the 36 pb _1 data sample acquired 

123 by the experiment in the 2010 [T5]. The Xc states are reconstructed through their 

124 radiative decay Xc J 1^1 with the J/ip decaying into two muons J/if> —> 

125 The results are shown in FigUJ compared with the CDF measurement and with two 

126 theoretical models, the color-singlet (blue) and NRQCD approach (red). 

127 6 Double charm production 

128 Recently, both the production of double quarkonium and also the associated pro- 

129 duction of quarkonium together with open charm have been suggested as probes of 

130 the production mechanism. In the pp collisions also other mechanisms, as the DPS 

131 [71 El E] , can be involved in the production and their contribution can be estimated 

132 with respect to the Single Parton Scattering (SPS) [UEEE]. 

133 Such a study has been performed at LHCb through the measurement of the double 
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Figure 7: Relative cross section cr(Xc) /cr(J '/?/>) in bins of <//■?/> transverse momentum. 
The results (black points) are compared with the CDF measurement [16] and with 
the color-singlet and color-octet prediction (respectively the blue and red area.) 



134 J/if>, taken from Ref.|17]. and J/ip production associated with an open charm hadron 

135 (such as D°, D + , D+ and A+) [18J, with the 2010 and 2011 datasets (respectively 37 

136 pb _1 and 355 pb _1 ). 

137 The double J/ip production cross section has been measured reconstructing the 

138 two J/tj) mesons in their decay to two muons. Both the J/ip mesons have been 

139 required to have rapidity and transverse momentum lying respectively in the ranges 
ho 2 < y < 4.5 and p T < 10 GeV/c. 

ui The signal yield is determined fitting the invariant mass distribution of the first muon 

142 pair in bins of the second muon pair and correcting the number of signal events by 

H3 the total efficiency. 

144 The total efficiency is factorized in three different terms 

J,ot sel&zreco&cacc C ^ID trg / R \ 

E j/i>j/i> — E J/ipj/i> x E j/tpj/i> x j/ipj/ipi \°) 

145 where £j e }^jf^ 0&cacc is the acceptance, selection and reconstruction efficiency, ^fRj,^ 
He is the efficiency of the muon identification and sj/^j/^ is the trigger efficiency. To 
H7 take into account the distortion due to the unknown J/ip polarization £ s j/ l ^jJ^ acc is 
us a function of the J/ip cos 8 where 9 is the angle between the // + in the J/ip center of 
H9 mass frame and the Lorentz boost from the laboratory frame to the J/ip frame. 

150 The corrected invariant mass distribution of the first muon pair in bins of the second 

151 muon pair is shown in FiglHl in three bins of J/ip transverse momentum in a particular 

152 bin of rapidity. 
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Figure 8: Invariant mass distribution of the first muon pair in bins of the second 
muon pair for the double J /if) production. 

153 The double J /if) cross section is estimated to be 

"W* = £ g2 (J/^V^-) = (5A ± 1 '° (Stat) ± 1A ( SyS ^ nb ' W 

154 where NJJ^j,^ is the efficiency corrected signal yield, C = 37 pb -1 is the integrated 

155 luminosity and B(J /if) — > /x + is the branching ratio of the J /if) decay into a muon 
we pair. 

157 The uncertainty are respectively statistical and systematic. The main contributions 

158 to the systematic uncertainty come from the tracking and trigger efficiency and from 

159 the unknown J /if) polarization. 

wo The experimental result obtained by LHCb has been compared with the theoretical 

lei contribution calculated in color-singlet model, from the SPS and the DPS. The two 

162 contributions, estimated in the LHCb acceptance, are listed together with the related 

«3 uncertainties, in TabJT] [17J. The sum of the two contributions is in agreement with 

164 the experimental value reported in EqJEJ although the uncertainties on the theoretical 

«5 expectations are too large to draw a definite conclusion on the production mechanism, 
we The production cross sections of a J /if) meson associated with an open charm 

167 hadron, D°, D + , or A+ have been measured using 355 pb _1 out of the 2011 

168 datasets. As control channels, the cc events with two open charm hadrons recon- 

169 structed in the LHCb acceptance have been also studied. The J/tf>, D°, D + , Df and 

170 hadrons have been reconstructed through the following decays: J /if) — > fi + [i~, 
in D° -»■ Tt+K-, D+ -»■ TT+7T+K-, Df tt+K+K-, A+ -»■ pn + K~. " 

172 In Figj9] (taken from Ref.[18j) the results for the production cross section of the J/ipD 
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Model 


Cross section (nb) 


Uncertainty 


Single Parton S. 


4.15 


30% 


Double Parton S. contribution 


2 


50% 



Table 1: SPS and DPS contribution to the double J/ip production cross section 
estimated in the LHCb acceptance. The theoretical uncertainties are also listed. 
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Figure 9: J/if> production cross section 
with associated open charm. The exper- Fig ure 10 ; 
imental results (black points) are com- °P en charm 
pared with the gluon-gluon fusion expec- °P en charm cross section ratio compared 
tation (yellow and green areas [I El 13). with the theoretical expectation com- 
puted with the DPS approach HH IS]- 



173 processes are shown, compared with the theoretical expectation estimated with gluon- 

174 gluon fusion model. In FigJTU] (from Ref.|18j) the ratios of the product of the prompt 

175 open charm cross sections and the double open charm cross section show a good 
ne agreement with the theoretical expectation from the DPS, assuming the effective 
177 cross section measured in multi-jet events at Tevatron [19]. In the plots the inner 
us error bars represent the statistical uncertainties while the outer error bars are the 
179 sum in quadrature of the statistical and systematic uncertainties. The main contri- 
lso butions to the systematics are coming from the trigger efficiency, the unknown J /if) 
i8i polarization and the luminosity. 



182 7 Conclusion 

183 LHCb has provided many contributions in understanding the quarkonium production 

184 mechanism in particular with the measurement of the J /if) and Xc production. Also 

185 the double J /if) and J /if) associated with an open charm hadron production has 
we been investigated, which provide a useful test for the SPS and DPS approach. The 
187 measurement of the polarization of the prompt J / if) component is ongoing and it will 
las provide a critical test for the validity of the color-singlet and color-octet models. 
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